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ACCRETION OUTBURSTS IN CIRCUMPLANETARY DISKS 
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Space Telescope Science Institute, Baltimore, MD 21218, USA 

ABSTRACT 

We describe a model for the long term evolution of a circumplanetary disk that 
is fed mass from a circumstellar disk and contains regions of low turbulence (dead 
zones). We show that such disks can be subject to accretion driven outbursts, 
analogous to outbursts previously modeled in the context of circumstellar disks to 
explain FU Ori phenomena. Circumplanetary disks around a proto- Jupiter can 
undergo outbursts for infall accretion rates onto the disks in the range M infa n « 
10 -9 to 10 MqVT -1 , typical of accretion rates in the T Tauri phase. During 
outbursts, the accretion rate and disk luminosity increases by several orders of 
magnitude. Most of the planet mass growth during planetary gas accretion may 
occur via disk outbursts involving gas that is considerably hotter than predicted 
by steady state models. For low infall accretion rates M infa n < 1CT 10 M yr _1 
that occur in late stages of disk accretion, disk outbursts are unlikely to occur, 
even if dead zones are present. Such conditions are favorable for the formation 
of icy satellites. 

Subject headings: accretion, accretion disks -- magnetohydrodynamics (MHD) 
- planets and satellites: formation -- protoplanetary disks -- stars: pre- main 
sequence 
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1. Introduction 

According to recent models, after a planet forms in a circumstellar disk and its mass 
reaches that of the orde r of Neptune's, tidal forces from the pla n et open a gap in the disk 
JLin fc Papaloizoul ll984J Il986t b'Angelo. Henning & Klevl I2OO2I : bate et all liooih . Mate- 
rial continues to flow from the circumstellar disk through the gap and onto the planet 



2003 



( Artymowicz fc Lubowll996 : Kley 1999 : Lubow. Seibert fc ArtymowiczlligggHD'Angelo. Henning fc Kiev 



Lubow fc D'AngelolbOOGuAyfiffe fc BateH2009l ). The planetary radius is much smaller 



than its Hill radius and, since gas flows into the Hill sphere with some angular momentum, 
a circumplanetary disk forms. If modeled as a standard viscous accretion disk, the circum- 
planetary disk is tidally truncated by the star at a radius of about 0.4 times the Hill radius 
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( IMartin fc Lubowll2011al ). Nearly all studies of the global structure of circumplanetary disks 
have assumed them to be fully turbulent. In this Letter, we explore the long-term conse- 
quences of dead zones in circumplanetary disks, regions that have no turbulence near the 
midplane. In particular, we show that these disks can undergo outbursts. 

The structure of a circumplanetary disk depends upon the rate of angular momentum 
transport across the disk. The angular momentum transport allows mass to move inwards 
and be accreted onto the central planet. Angular momentum transport is driven by tur- 
bulence, of which there are two main sources. First, the magneto rotational instability 
(MRI ) drives turbulence if th e gas is sufficiently well ionized to be coupled to a magnetic 
field (IBalbus fc Hawleylll99ll ). In the hot inner parts of the disk, the field is well coupled 
due to thermal ionization and the MRI drives turbulence. However, further out where the 
disk is cooler, the low ionization fraction causes magnetic field instabilities to be suppressed 
( iGammid 1 19961 ; iGammie fc Menoul Il998l ; iTurner fc Sand 120081 ) . The disk is said to form a 
dead zone at the midplane because it is not turbulent in this region. However, the disk 
becomes layered because external radiation suc h as cosmic rays or X-rays can i onize the sur- 
face layers, allowing them to be turbulent (e.g.. lGlassgold. Najita fc Igeall2004l ). The second 
type of turbulence is driven by gravi tational instability if the disk acquires enough mass 
rtPaczvnskilll~978l ; lLodato fc Ricel 120041 ) . 



Young circumstellar disks are thought to undergo out bursts (FU Ori phenomena) as a 



consequence of a sudden change in the level of turbulence (JGammid Il996t lArmitage at al. 



200ll ; IZhu et al.l 120091 ; iMartin &: Lubowl l2011bf ) . Because a dead zone inhibits inward flow 
onto the star, mass should build up within the circumstellar disk. With sufficient mass, the 
dead zone becomes gravitationally unstable. The instability causes turbulent heating and 
increased ionization of the gas, which can trigger the MRI. This triggering causes a sudden 
increase in the disk turbulence, resulting in an accretion outburst. After the outburst, the 
remaining disk cools sufficiently for a dead zone to form again and the outburst cycle repeats 
period ically. This process, the g ravo-magneto (GM) instability, can be understood as a limit 
cycle ( IMartin &: Lubowl 1201 lbl ). The outbursts can be explained as transitions between 
steady state disk solutions in a diagram that plots the disk effective temperature versus the 
total surface density at a fixed radius. This graph is analogous to the well-known S-curve 
used to explain dwarf-nova outbursts. There is a range of accretion rates for which no steady 
disk solution exists and for these rates the disk is GM unstable. 



Similar outbursts could also occur in circumplanetary disks. IMartin fc Lubowl ( 12011al ) 
suggested that dead zone formation is even more favorable in circumplanetary disks than 
in circumstellar disks because the surface densities may be higher, while the temperatures 
remain low enough to avoid thermal ionization. In this work, for the first time we model 
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the global structure of a layered circumplanetary disk over long timescales for a variety of 
parameters. In Section [2] we describe the disk model and the disk evolution. In Section [3] 
we discuss the observability of a circumplanetary disk. Finally, in Section H] we discuss 
implications of the layered circumplanetary disk model. 



2. Circumplanetary Disk Models 



2.1. Model Description 



We apply the l ayered disk mod el in itially described in lArmitage at al. I (1200 ll ) and 
further developed by lZhu et al.l (120091 ) and lMartin fc Lubowl (J2011bl ). Material falls onto the 
disk at a constant rate Mmfaii and is accreted onto the central planet through the inner disk 
boundary at a rate M. The model accounts for the effects of disk turbulence by a viscosity. 
The disk vertical structure consists of layers: outer magnetically turbulent surface layers 
on the upper and lower disk surfaces and a midplane layer that can be either nonturbulent 
(dead) or turbulent (active) due to the effects of disk magnetic fields and/or self-gravity 
The disk temperatures in these layers are accounted for by a T — r relation that is based 
on dust opacity. Disk heating is taken to be only due to accretion. External heating by 
the star or planet is ignored. Heating by the star may be ignored for planet orbital radii 
> 1 AU, since its contribution to the disk temperature is much less than that required for 
thermal ionization ~ 800 K. The disk is taken to be in Keplerian orbit about the planet. 
The model does not account for the effects of pressure or self-gravity on the disk rotation. 
The planet is assumed to accept whatever gas is accreted through the disk inner boundary 
without significant ex pansion into the disk, as in the run-aw ay phase of gas accretion of the 
core accretion model (ID'Angelo. Durisen. fc Lissauer Il2010h . We take the mass of the planet 
to be fixed and do not account for its mass increase through gas accretion. 

The disk has a total surface density E(i?, t), midplane temperature T c (R,t) and surface 
temperature T e (R,t) for radius R and time t. The MRI turbulent surface layers have total 
surface density T, m (R,t), temperatu re T m (R,t), and viscosity u m (R,t) that is parameterized 
with the IShakura &: Sunyaevl ( 119731 ) a m parameter. The external radiation is assumed to be 
able to penetrate to a constant surface density E crit in total for the disk upper and lower 
surface regions. If the surface density is greater than this critical value, E > E crit , then 
E m = Ecrit- For smaller surface density, E < E cri t, the disk is MRI active at all heights, so 
that E m = E. 

A complementary midplane layer exists if the surface density is larger than the critical 
value, E > E crit . It has surface density E g = E — E crit and midplane temperature T c . This 



-4- 



layer is MRI active if the midplane temperature is greater than the critical, T c > T cr i t , where 
T crit is the temperature for sufficient thermal ionization. However, for lower temperatures, 
T c < Tent, a dead zone forms. There is then no turbulence near the midplane unless it is 
massive enough to be self-gravitati ng. The condition for self-gravity is taken to be that the 
Toomre parameter (jToomre Ill964l ) is smaller than its critical value, Q < Q C mt = 2. In this 
case, a second viscosity term is included in the complementary layer. 

We consider the case of a Jupiter-like planet of mass, M p = lMj, orbiting a solar 
mass star, M = 1M , at a distance of a = 5.2 AU. We solve the accretio n disk equations 



numerically on a fixed me sh that is uniform in Ri with 120 grid points (e.g.. lArmitage at al. 



200ll ; iMartin et al.l 120071 ). The inner boundary has a zero torque condition at the radius of 
Jupiter, R = IRj, so the mass falls freely onto the planet. At the outer boundary we 
approximate the tidal torque from the star wi th a zero radial v e locity boundary condition 
at -R out = 0.4 i?n, equivalent to equation 70 of IMartin fc Lubowl (j2011al ). Circumstellar gas 
accretes onto the circumplanetary disk at a constant rate Mi n f a n at a radius -R a dd- 

The specific angular momentum of the accreting gas entering the planet's Hill sphere 
from the circumstellar disk has yet to be determined. One estimate is obtained by assuming 
that the gas enters near the Lagrange points, L% and L2, and is nearly at rest relative to 
these points as it enters. By conservation of angular mo mentum, the circularization radius 
for infalling material is estimated to be about 0.33 i?n (IQuilleru Il998t Estrada et al.l 120091 ) 
and we take this radius to be i? a dd = 0.33 Rh- We have also considered cases of i? a dd = 
0.1 Rh and found that the results are about the same. This outcome is consistent with the 
results for a fully turbulent disk by IMartin fc Lubowl (j2011al ) that show the circumplanetary 
disk structure is not very sensitive to the angular momentum distribution of the incoming 
circumstellar gas. Initially we take the surface density to be very small everywhere and allow 
material to build up. We run the simulation until either the disk reaches a steady state or 
the mass accretion onto the central planet shows an outburst pattern that repeats itself. 

The properties of the turbulent surface layer have not been well determined. There have 
been recent attempts to determine its depth, S cri t, by calculating in detail the ionization 
balance of the disk due to external sources of radiation and subject to various processes, 
such as ambipol ar diffusion and the pr e sence of PAH and dust that tend to supp ress the 



ionization (e.g., iBai fc Goodman 1 120091 ; IPerez-Becker fc Chiang 1 12011c iBai 1120111 ). When 



applied to protostellar disks, these calculations result in accretion rates that are lower than 
observed in typical T Tauri stars. The discrepancy is currently unresolved. Rather than 
attempt to fold in such an approach in this Letter, we follow Armitage et al. (2002) and 
Zhu et al. (2009, 2010) and regard T, c r \t as a free parameter. Since ty pical T Tauri accretion 
rates suggest that S cri t > 10 gem -2 (IPerez-Becker fc Chiang 1 1201 ll ). we consider values in 
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the range 100 g cm 2 > S cr i t > 10 g cm 



Following Armitage et al. (2001), we take T crit 
viscosity a m parameter is found to be ~ 0.01 — 0.4 fe.g.lHartmann et al. 



800 K. Various arguments suggest the 



1998: 



Fromang et al. 



2007l : lGuan et al.lbood ; E>avis et alJl2009Hzhu et al.ll2007l : lKing. Pringle fc Livioll2007h . With 
the uncertainty in the value, we consider cases with a m = 0.01 and 0.1. 

During times of planet formation, the accretion rate onto the circumplanetary disk, 

Mnfaii, is of order the overall circumstellar disk accretion rate (Bate et al. 2003; Lubow fc 

D'Angelo 2006; Ayliffe & Bat e 2009) that we take to be of order 10~ 8 M yr" 1 JValenti. Basri fc Johns 
19931 ; iHartmann et al.lll998l ). typical for the T Tauri phase. In this Letter, we consider a 
range of values for a m , M infall , and E crit . 



2.2. Model Results 

In Table [1] we show the outcome for a set of model parameters. If the model does 
not reach a steady state and outbursts occur, then we tabulate the time interval between 
outbursts, ti n t, and the duration of outbursts, £ out . For a m in the range 0.01 — 0.1 with 
S C rit ;$ 100 gem -2 , outbursts are possible for accretion rates onto the disk of M infal i = 10~ 9 
to lO _7 M yr _1 . They repeat on a typical timescale of 10 4 to 10 5 yr and last for a few to 
several years. We caution, however, that in the outburst state the disk thickness is formally 
of order the disk radius. Consequently, the idealization of a disk geometry in Keplerian 
rotation at this stage is questionable. In the quiescent state, the disk thickness is a few 
tenths of the disk radius, as ex pected by equat i on 4 of iMartin fc Lubowl (|2011aJ) and the 
three-dimensional simulations of lAyliffe fc Batd (120091 ). We find that the disk mass reaches 
values that are not very small compared to the planet mass. For example, model R3 reaches 
values as high as about 0.3 Mj. The Keplerian disk assumption is then not well satisfied at 
all times. Improved mult i- dimensional modeling should be carried out in future work. 



M max , of about 10" 4 M 



q yr x and a 



The results indicate a maximum accretion rate 
quiescent accretion rate of about lO _9 M yr _1 . Consequently, accretional heating of the 
gas increases substantially during an outburst. The mass accreted during an outburst is 
somewhat greater than the mass accreted between outbursts, typically 30% to 40% greater. 
The disk is stabilized against outbursts with a larger active outer layer surface density, E crit , 
or stronger turbulence a m . For smaller values of S crit , outbursts become more problematic. 
Model R7 with S cr j t = 10gcm~ 2 has a longer interval between outbursts and involves an 
accreted mass during outburst that is comparable to that of the planet which suggests a 
more detailed treatment is required to determine whether outbursts occur in this regime. 



-6- 



The detailed ionization models discussed in Section [27X1 obtain even lower values of S crit , for 
which our models are subject to this limitation. 

The instability cycle can be understood through a state transition diagram that plots 
effective surface temperature (or equivalent stea dy state accretion rate) a s a function of disk 



surface density at some fixed radius in the disk (IMartin fc Lubowll2011bl ). In Fig. [T] we show 
the limit cycle for model R3 in the such a diagram at a radius R = 0.07 Rn. The steady 
ac cretion rate, M s , correspo nds to the surface temperature of the disk given by equation 18 



in 



Martin fc Lubowl (l2011bl ). The thick solid lines show the two types of analytical steady 
state solutions. The first is fully MRI active and this has two branches, one with T > T cr ; t 
(upper branch) and a second with S < S crit (lower branch). The second type of steady 
state is the gravo-magneto (GM) solution that has a self gravitating dead zone and MRI 
active surface layers. The dotted line shows the nonsteady dead zone branch that connects 
the lower MRI branch to the GM branch. The numerical simulations reveal a cycle on this 
S — M s diagram. The simulation results plotted in Fig. [1] generally track along the dead 
zone and steady state branches. The simulation track moves to the right from the dead zone 
branch to the GM branch. This branch is plotted as a short line segment at the right end 
of the horizontal line. After outburst, the simulation backtracks somewhat to the left and 
slightly upward, and then advances directly upward to the upper MRI branch. This upward 
transition marks the onset of the outburst. 

Between the dead zone/GM branches and the upper MRI branch there is a range of 
accretion rates M s for which no steady state is possible and an outburst results. This unstable 
regime is plotted as the shaded region. Model R3 lies in that range at the selected radius, 
since it has an infall accretion rate of 10 -8 M yr" 1 . The global range of unstable accretion 
rates covers a somewhat wider range than is shown on the plot, since other unstable rates 
may occur at other radii. 

The results in Table [1] are consistent with there being a limited range of infall accretion 
rates that result in outbursts. Model RIO undergoes outbursts while models R9 and Rll, 
having lower and higher infall accretion rates respectively, do not undergo outbursts, with 
all other input parameters being the same. 



3. Observability of a Circumplanetary Disk 

The average accretion rate onto the planet is equal to the rate of infall accretion onto 
the disk. Because the accretion through the disk is time dependent, so is the luminosity. 
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Fig. 1. — Limit cycle for model R3 at a radius of R = 0.07 i?n in the £ — M s state or E — T c 
diagram. T e is the disk surface temperature that is less than T c , the temperature at the 
disk midplane that controls the onset of MRI via thermal ionization. The fully turbulent 
(no dead zone) solutions involving MRI have two branches. They are plotted as the sloped 
heavy line with a gap. The lower branch (below the gap) is ionized throughout by external 
sources of radiation and the upper branch is thermally ionized. The steady gravo-magneto 
(GM) solution is shown as a short segment on the right hand end of the horizontal line. The 
dotted line shows where the nonsteady dead zone branch resides that connects the lower 
MRI branch to the GM branch. The thin line shows the simulated evolution in time that 
evolves in a counter-clockwise direction. It extends horizontally to the GM branch. The 
shaded region shows the range of accretion rates for which the disk has no steady solution 
and would trigger an outburst. 
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Table 1: The parameters of each model and the timescales of the outbursts if they occur. 
Column 2 is the a m parameter in the active outer layers, column 3 is the accretion rate onto 
the disk in units of M yr -1 , column 4 is the total surface density of the active outer layers 
(in gem -2 ), column 5 is the time interval between outbursts, column 6 is the duration of 
the outburst, column 7 is the mass accreted onto the planet during the outburst, column 8 
is the accretion rate during quiesence and column 9 is the maximum accretion rate onto the 
planet achieved during the outburst. 
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The luminosity of an accretion disk is 

GMM 



2R h 



(jPringld 119811 ) , where R m is the inner edge of the disk that is assumed to coincide with 
the radius of the central star or planet. For example, a disk around a solar mass and 
radius star with an accretion rate of 10~ 8 M yr -1 has an accretion luminosity of 0.15 L Q (in 
addition to the luminosity due to heating from the star). If a disk surrounds a Jupiter mass 
and (constant) Jupiter radius planet that has opened a gap in the circumstellar disk, the 
quiescent luminosity of the circumplanetary disk is 0.0015 L or only 1 percent of accretion 
luminosity of the circumstellar disk. 

If the disk is observed during an outburst with a high accretion rate of lO _4 M yr _1 , 
then the luminosity is ~ 15 L Q (see Tabled]). However, the disk will be in outburst only for 
0.01% — 0.1% of time. It is possible that with the increased luminosity even for short periods 
of time, forming planets may be more easily detected. Synoptic surveys may be capable of 
detecting such events in the future. 



4. Discussion and Conclusions 

Circumplanetary disks can contain regions with low turbulence (dead zones) because 
they have high enough surface densities to prevent ionization of the disk midplane by external 
sources of radiation, but remain cool enough to avoid thermal ionization. As a result of dead 
zone formation, the disks are unstable to the gravo-magneto (GM) instability in which a 
sudden turn-on of the MRI triggers a rapid accretion event (see Fig. 1). We have found that 
such outbursts are possible for infall accretion rates onto the disk in the range of 10 -9 to 
10~ 7 M yr -1 , typical of accretion rates in the T Tauri phase of young stars, when gas giant 
planet formation may occur. The quiescent accretion rates onto planets are ~ 10~ 9 M Q yr -1 , 
but during outbursts the rates increase up to ~ 10~ 4 M yr -1 , making them more observable. 
The outbursts typically occur every ~ 10 4 to 10 5 yr and have a duration of only several years. 
Somewhat more gas is accreted during outbursts than between them. 

The results have implications for the nature of gas accretion onto a planet during the 
run-away gas accretion phase of the core accretion model (see D'Angelo et al. 2010). The 
higher entropy we expect to be associated with the gas accreted onto a planet in outbursts 
may have implications on the initial state of a newly formed planet. The higher entropy may 
favor more luminous "hot start" planet models (c.f. Marley et al 2007). 

There are several approximations in the layered disk model we have considered that 
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should be improved upon in future work. Some of these were discussed in lMartin fc Lubow 



( J2011bl ) in the context of circumstellar disks. The present treatment is based on a one- 
dimensional model of a disk in Keplerian rotation. As noted in Section 12.2^ in the case of 
circumplanetary disks, the disk thicknesses and masses can result in structures that may 
not be well approximated by a Keplerian disk. Multi-dimensional simulations that include 
pressure and self-gravity would be of benefit in understanding the evolution of the outburst 
and the accretion onto the planet. 

Lower accretion rates (< 1O~ 1O M yr _1 ) provide lower disk tempe ratures that would be 



compatible with the survival of icy satellites at thei r current locations (jCanup fc Wardll2002 



Mosqueira fc Estrada 1 120031 ; IWard fc Canupll2010l ). Such rates can occur at late stages of 



disk evolution. The Jovian circumplanetary disk could have contained a dead zone in the 
regular satellite region, provided that external sources of radiation penetrated a disk surface 
density of S crit < 10 gem -2 . It is unlikely that outbursts will occur at this stage even in 
the presence of a dead zone, since the time required for a disk to reach a self-gravitating 
state is longer than the disk lifetime. Therefore, although outbursts could adversely affect 
the survival of icy satellites, they are unlikely to occur at these late stages. 
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